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Expressions are developed which permit calculation of the temperature-dependent 
thermal conductivity of a cylindrical specimen from the experimentally determined tem- 
perature distributions on its surfaces and from the heal How through a central circular area 
at one end of the cylinder. Numerical factors are tabulated which permit rapid hand 
computation of thermal conductivity values from experimental data for the case of nxially 
symmetric parabolic boundary conditions on the flat surfaces and a linear longitudinal 
temperature distribution on the convex surface. Applications to several published methods 
of thermal conductivity determinations are shown and examples are given. 



1. Introduction 

Tn a recent paper by Glaser el al. [1], ] a mot hod 
of measuring the t hernial conductivity of solids was 
described in which the principal heal flow was 
parallel to the axis of a small circular disk-shaped 
specimen, hut in which significant radial heat 
flows also occurred. Assuming axial symmetry, the 
authors used a least-mean-squares parabola to 
represent the measured radial temperature profile 
on each face of the specimen and assumed a linear 
longitudinal temperature distribution along the 
cylindrical surface of the specimen. Utilizing these 
boundary conditions, longitudinal temperat ure gradi- 
ents at the cooler face of the specimen were deter- 
mined by numerical solution of a finite difference 
equation using a digital computer. The thermal 
conductivity of the specimen was calculated from 
the average computed temperature gradient at the 
cooler surface and the measured heat flow leaving 
this surface. 

A somewhat similar heat flow analysis is being 
used in the Heat Transfer Section of the National 
Bureau of Standards in conjunction with thermal 
conductivity measurements employing a guarded 
steam calorimeter for the measurement of heat 
flow through the central 3-in.-diam area of a 6-in.- 
diam disk-shaped specimen [2]. The method is one 
of chiefly longitudinal heat flow, but corrections are 
occasionally necessary to account for radial heat 
flows. In order to effect such corrections, an exact 
analytical solution was obtained to the boundary 
value problem in which parabolic temperature dis- 
tributions having axial symmetry are assumed on 
the two plane end surfaces, and a linear longitudinal 
temperature distribution is assumed on the convex 
surface, of a right circular cylinder. 

Subsequent to the development of this analysis 
at NBS, a paper was published by Hoch et al. [3], 
in which a somewhat similar, although less general, 



1 Figures in brackets indicate the literature references at t he end of t his paper. 



analytical solution was derived. In the latter paper, 
the heat conduction equation was solved for the 
case of a finite cylinder with the convex surface at 
a uniform temperature and the same axially sym- 
metric parabolic temperature distribution on the 
two laces of the cylinder. The thermal conductivity 
of the specimen was calculated from the computed 
values of the temperature gradient at the center of 
the faces of the cylinder and of the heat (lux radiated 
from these surfaces at the center, as determined by 
their temperature and emissivity. 

In view of the general interest in this heat flow 
problem, it was decided to present, prior to formal 
publication of the paper [2] describing the NBS steam 
calorimeter apparatus, the mathematical analysis de- 
veloped for it. The analysis given first treats the 
case of a right circular cylinder with arbitrary tem- 
perature boundary conditions, and then develops 
solutions for specific boundary conditions applicable 
to the determination of thermal conductivity. 

2. Mathematical Development 

Consider a homogeneous, isotropic, opaque, solid 
right cylinder of radius, />, and thickness, /. In 
general, the thermal conductivity of the cylinder 
material may vary with temperature. If so, the 
steady-state heat flow equation in cylindrical co- 
ordinates, assuming axial symmetry, is 



V-(kVv)-- 



1 d_ 

r dr 



k(v)r 






= 0, (1) 



where temperature above an arbitrary datum is 
denoted by the symbol v, the temperature-dependent 
thermal conductivity by k(r), and the radial and 
longitudinal coordinates by r and z, respectively. 
Equation (1) can be reduced to a simpler form by 
introduction of a new variable, y, defined by the 
relation 



y= 



F.(>'^ 



(2) 



215 



where k* is the value of k(v) at v = 0. 2 
(1) reduces to [4, 5] 



dr 2 



dr ' d2 



Utilizing (2), 



(3) 



The variable, ?/, is essentially a potential whose 
gradient is proportional to the heat flux [6]. 



2.1. General Boundary Conditions 
Assume the following boundary conditions: 
2 = 0<r<b y=j\r), 

2=1 o<r<b y=g(r), 



(4a) 
(4b) 



0<z<7 r=b 



y=J(b) + \g(b)-J(b)] 



I 



+h(z), (4c) 



where /(r) and g{r) are arbitrary potential distribu- 
tions on the flat surfaces of the disk, and where, in 
order to insure continuity, we shall require that the 
function h (z) vanishes at the flat surfaces of the 
disk, so that h(0)=h(l)=0. The solution is [7]: 3 

y=f(b)+[g(b)-J(b)] f+ge/o (<*« 



X 



1 A n sinh [a n | \+B n sinh ( a n —^-) V 



+J2CJ Q 



kirr \ 
~T) 



sin 



Kr)' 



(5) 



where J and I Q are, respectively, the ordinary and 
the modified Bessel functions of the first kind; the 
a n 's are the positive roots of J (a n )=0. In order to 
satisfy the boundary conditions, the coefficients 
must be: 



A 



¥J\(a n ) sinh (aj/b) 



I { 9(r)—g(b) }rj (a n r/b)dr, 



B n - 



b 2 J\{a n ) sinh (a n l/b) 



j> 



and 



(6) 
(7) 



Cr 



Uc 



^J> ( 



(¥) 



z) sin f — j-\ dz. 



(8) 



Equation (5) can be used to calculate the steady- 
state potential distribution in a finite cylinder corre- 



2 This is a special case of the more general relation dy=Ck(v)dv, where C is an 
arbitrary constant. 

3 This specific problem is not handled by Carslaw and Jaeger; the form of solu- 
tion and the necessary integrals are given in the chapter cited. 



sponding to arbitrary potential boundary conditions, 
g(r),f(r), and h (z). If the indicated integrations of 
(6), (7), and (8) cannot be obtained analytically, 
values can be determined by numerical integration. 

2.2. Parabolic Boundary Conditions 

The boundary conditions used by Glaser et al. [1], 
Robinson et al. [2], and Hoch et al. [3] are all of the 
form 

,2 



f(r)=Y +E 01 - 2 , 



op 



and 



g(r)=Y 1 +E 1 ^ 2 > 

h(z) = 0, 



(9) 



where Y and Y\ are the potentials (corresponding 
to y) at the centers of the faces of the disk, and E 
and Ei are the potential differences between the edge 
and the center of the disk at z = and z = l, respec- 
tively. With these boundary conditions, (5) becomes 

y=(Y +E )+(Ti-Y +E 1 -E ) " 



I 

J {a a r/b) 



X 



{ 



jfrio&TiCa,,) sinh (aj/b) 
1 sinh (a n £)+E sinh (a n ^) | ■ (10) 



3. Calculation of Thermal Conductivity 

The thermal conductivity of the disk is to be 
determined from the measured heat flow leaving a 
central portion of the disk of radius, r=a, and from 
the observed potential distributions on the surfaces 
of the disk. A mathematical expression giving the 
appropriate relationship between these quantities is 
derived by determining the total heat flow through a 
circle of radius a, at z=0: 

e-*rj>*(s) ,.*-«•// (S),.*- 

(11) 

Performing the indicated differentiation and integra- 
tion on (5), 

x{A-5 B cosh(«4)}+?§C7^(^)] 

(12) 

Equation (12) can be rewritten in terms of the 
parabolic boundary conditions (9) discussed in sec- 
tion 2.2. If this is done, the heat flow through a 
circle of radius, a, becomes 
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in which 



Q= Ta ' ±- ^-Yo+E&i-E&o] 



*.=i-i6 ± ,^ a ;i h ).., 



(l/b) 



and 



n=i(a/b)aUi(a, l ) sinh (aj/b) 
{lib) 



%=l-16± r ± (a fJ h [, X 



n=i (alb)alJi(a,,) Isttlh (aj/b) 



(13) 



(14) 



(15) 



are factors which depend only on the geometry of t lie- 
disk and on the fraction of the disk from which the 
heat flow is measured. If all of the other quantities 



are known, eq (13) can be used to compute the 
thermal conductivity, k*. In order to compute the 
values of the potential, y, for use in eq (13), some 
assumption must be made regarding the temperature- 
dependence of the thermal conductivity of the disk 
(see eq (2)). This is discussed in section 4. 

In the special case, l/b = 0, (14) and (15) reduce to 

n =i (a b)otiJ \(ot n ) 2 \o/ 



In the special case, a/b=0, used by Hoch et al. [3], 
Table 1. The coefficient ^i as a function of l/b and alb for parabolic boundary conditions 
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. 531 186 
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. 52424 
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1.0 
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.81414 
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. 82340 
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. 87724 
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2.0 
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. 9(5017 
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2. 1 


2.6 
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. 98526 


.'.'MilS 
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3.0 
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4.0 
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and 



* -1-8 s _ 1__ X __JM___ 

£=i a*Ji(aJ sinh (a n Z/6) 



Vi-8S^ X; m 



„=i aJJ^a,,) tanh (aJ/6) 
For large values of l/b, 



and 
where 



Z=i (a/b)a 3 n Ji(a n ) 



(17) 

(18) 

(19) 

(20) 

(21) 



is a function of a/b only. 

Using a high-speed digital computer, numerical 
values for ^ and % have been calculated for a range 
of values of a/b and l/b and are given in tables 1 and 2. 
Representative values are plotted in figure 1. The 
series for ^ is rapidly convergent, due to the quan- 
tity, sinh (a n l/b), which appears in the denominator. 
The series for ty , however, converges rather slowly 
and it was necessary to sum at least the first one 
thousand terms of this series in order to obtain the 
number of significant figures aiven in table 2. 



I - 



a/b = I.O ^- 



m 



W 



*--- a/b = 0.5^/ 




Figure 1. The coefficients ^ and ^i for the case of parabolic 
boundary conditions; solid lines, ^o; dashed lines, ^i. 



For vaJues of l/b larger than those given in tables 1 
and 2, ^ is substantially unity and ^ may be com- 
puted using (20) and the values of <£ given in table 3. 

Table 3. The function 3> which is used for computing values 
of ^o corresponding to large l/b 
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0.3 


2. 07864 


0.4 


2. 03028 


0.5 


1. 99630 


0.6 


1. 88502 


0.7 


1. 78393 


0.8 


1. 65885 


0.9 


1. 50226 


1.0 


1. 29410 



4. Effective Mean Temperature 

The thermal conductivity of the disk at v=0 can 
be computed from (12) for the general case or from 
(13) for the case of parabolic boundary conditions. 
The variable, v, is the temperature measured from 
an arbitrary reference temperature, say T*. Thus 
the conductivity obtained, k*, corresponds to the 
reference temperature J 7 *. The variable, y, as men- 
tioned previously, depends on the temperature- 
dependence of the thermal conductivity of the disk. 
Two particular cases of interest are discussed here. 

4.1. Thermal Conductivity Having Linear 
Temperature-Dependence 

For many materials the thermal conductivity can 
be assumed to vary linearly with temperature over a 
moderate temperature range, such that 



then, from (2), 



k(v)=k*[l+Pv]) 

I & 2 



(22) 
(23) 



The boundary conditions must be written in terms 
of y and then substituted into the appropriate 
equations. 

For the case of parabolic boundary conditions, the 
following relations hold: 



Y =V +% VI 



E =D +^ (Dl+2D V ) 



Y l =V 1 +£v l 



(24) 



(25) 



where F„= T - T* and V^T.-T*, T and T x being 
the temperatures at the center of the disk at 2 = 
and z = l, respectively, and Z? and Di are the tem- 
perature differences between the edge and the center 
of the disk at z=0 and z=l, respectively. 
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In order to compute the thermal conductivity of 
the disk, eqs (24) and (25) must be substituted into 
(13). The resultant expression, however, involves 
/3, the temperature coefficient of the thermal con- 
ductivity, a quantity which may not he known. 
This dependence can be eliminated by proper choice 
of the reference temperature, T*. II* T* is defined as 



2[T 1 -T +D{* 1 -/) l) * l) \ 



(26) 



then eq (13) reduces to 

Q=ira 2 j [Ti-To+D^-D^o], (27) 

which is the expression that would be obtained for 
the case of a specimen having constant thermal 
conductivity, k*. For the special case in which 
there is no radial heat flow, Di=D =0, and (26) 
reduces to T*=(2 T 1 +T )/2 [8]. 

4.2. Thermal Resistivity Having Linear 
Temperature-Dependence 

For many crystalline dielectric solids, the thermal 
resistivity, i.e., the reciprocal of thermal conduc- 
tivity, varies substantially linearly with temperature 
over a considerable temperature range. In this 
case, the thermal conductivity can be written as 



then, from (2), 



l+yv 



y=- In (I+70). 

7 



(28) 



(29) 



The relation between y and o can be better seen by 
expanding (29) in series form, 



y=v (l 



yv | (t?0 2 

"2"*" 3 



...> 



(m 



provided |y#|<^L 



In the present case the coefficient y is not so 
amenable to elimination by choice of a particular 
reference temperature as was the coefficient fi in 
section 4.1. In computing thermal conductivity 
values for a material whose thermal conductivity is 
expected to have a temperature dependence such as 
that in (28), one may first assume that \yv\ is suffi- 
ciently small to permit terminating eq (30) after 
the second term of the series; i.e., y=v—yv 2 /2. 
Comparison of this expression with (23) shows that 
the assumption that \yv\ is small is equivalent to 
assuming a linear temperature dependence for the 
thermal conductivity. Thus one may use (26) and 
(27) to compute a mean temperature and the ap- 
proximate thermal conductivity at that temperature. 
The first approximations to the thermal conductivity 
thus obtained for several tests at different mean 
temperatures can be used to compute values of the 



coefficient y which then can be used in (29) or (30) 
to obtain values of ?/, which can in turn be used in 
(13) to obtain more refined values for thermal 
conductivity. The resultant value of k* can be 
used to obtain better values of y, etc. The con- 
vergence is quite rapid in most practical cases, so 
that only one or two iterations will normally be 
required. 

5. Examples 

In order to illustrate the use of the mathematical 
developments presented in this paper, we now shall 
cite specific examples in which the factors % and 
^o given in section 3 for parabolic boundary con- 
ditions are utilized to compute thermal conductivity 
values from measured temperatures and a measured 
heat flow. One example will be given for each of the 
three methods mentioned in section 1 [1, 2, 3]. 

5.1. Method of Glaser et al. 

In the method described by Glaser et al. [1], the 
entire heat flow leaving the colder flat surface of the 
disk-shaped specimen is determined, i.e., in the 
nomenclature of the present paper, a/b = 1.0. We 
postulate the following data as having been acquired 
experimentally: 4 

a = b = 0A76 cm Z=0.16() cm 

#=9.60 W 
T =1319 °C A)=-218 deg 
7; = 1536 °C A= -338 deg 

Interpolation of table 1 and 2 for the case a/b =1.0 
and l/b = 0.336 yields, 



^=0.559 



*o = 0.372. 



Assuming linear temperature dependence, the 
thermal conductivity is computed using (27). Sub- 
stituting the values given above for the various 
parameters and solving for thermal conductivity, 

£* = 0.0198 W/cm deg. 

The mean temperature is computed from (26) as: 

r*=1370 °c. 

It should be remembered that, for a given specimen 
size and shape, the factors ^i and ^ remain constant. 
The only further information needed to compute 
thermal conductivity values for a given specimen 
is the heat flow and four temperatures (in effect). 
Thus the procedure described by Glaser et al. [1], in 
which each thermal conductivity value was obtained 
by numerical solution of a finite difference equation 
on a high-speed digital computer, can now be 
replaced by a fairly simple hand calculation using 
the values of % and ^ tabulated in this paper. 



4 This is an actual set of data which was made available through the courtesy 
of P. E. Glaser of Arthur 1). Little, Inc. 
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5.2. Method of Robinson et al. 

In the NBS steam calorimeter method [2], the 
total heat flow through only a central section of the 
specimen is measured. Let us assume that the 
thermal resistivity, rather than the thermal con- 
ductivity, of the specimen exhibits a linear temperature 
dependence; i.e., the thermal conductivity is of the 
form given in (28). We postulate the following 
data as having been acquired experimentally: 



a = 3.81 cm 6 = 7.62 cm 


£=2.29 cm 


#=542 W 




T =540 °C Z?o=- 


-20 deg 


Ti = 980 °C A = - 


-40 deg 



For the case a/b = 0.5 and l/b = 0.3, table 1 and 2 yield 

^ = 0.18481 ^o = 0.00520. 

Substitution into (26) and (27) yields 

&* = 0.0629 W/cm deg at T'* = 756.5 °C 

as a first approximation to the thermal conductivity 
at this temperature. If similar sets of data at other 
temperature levels are analyzed in this manner, the 
resultant thermal conductivity values can be used 
to compute values for y, the thermal resistivity 
temperature coefficient. Let us assume that this has 
been done and that a value of 7 = 9.0 X10 -4 deg -1 at 
757 °C was obtained. Using T* = 757 °C as a 
reference temperature, 



F =-217 deg 



Fi = 223 deg. 



Substitution of these values for V and V\ and the 
above-mentioned values for D Q , D u and y into (29) 
yields 

F = -241.4 deg £ T o=-25.1 deg 
Fi= 203.2 deg E x = — 33.8 deg 

Substitution of these values of Y , Y h E , and E x 
into (13) yields 

fc*=0.0621 W/cm deg at T* = 757 °C 

as the improved value for the thermal conductivity. 
If thought necessary, the improved values for the 
thermal conductivity can be used to compute an 
improved value for y for further refinement of the 
final thermal conductivity values. 

5.3. Method of Hoch et al. 

In the method described by Hoch et al. [3], the 
heat flux from the center of the flat surfaces of the 
disk is calculated from the temperature at that 
point and the total hemispherical emittance of the 
specimen, using the Stefan-Boltzmann radiation law. 
Since Hoch et al. derived an analytical solution for 
their case which is quite similar to that developed in 



the present paper for more general conditions, it is 
of interest to compare these solutions. Hoch et al. 
define a characteristic constant, K , which is given 
by 5 

K =2 £ tMhy/26) , (31) 

n =i a n Ji(a n ) 

where the notation has been changed to conform with 
that of the present paper. 

With the simplified boundary conditions assumed 
by Hoch et al. (Y l =T 0i E 1 =E ), the heat flux at 
the center of the surface of the disk is given by 



Kf 



) 



/ 



■E [% — *o 



(32) 



This expression is analogous to (13), except that it 
refers to the heat flux at a point rather than to the 
total heat flow through a circular area. For the 
case a/b = 0, with parabolic boundary conditions, the 
factors % and ^ are given by (17) and (18); the 
quantity [S^ — %] becomes 



[*i-*o]=8 S 2 w v r . ) 
n=i<x n Ji(a„) \_U\nh [a n 



lib) 



1 



which reduces to 



-Ml/b) S 



sinh (a n l/b)j 
tanh (a n l/2 b) 



Comparison of (31) and (34) yields 



(33) 



(34) 



(35) 



Since Hoch et al. included a numerical example in 
their paper, it is needless to present one here. Values 
for Kq can be computed using eq (35) and the appro- 
priate values of % and ^ from the first column in 
tables 1 and 2. Hoch and Nitti [9] more recently 
have published an expanded version of their previous 
paper [3] in which they tabulate values of K for 
several specimen shapes. 

6. Discussion 

For the case of parabolic radial potential distri- 
butions on the flat surfaces of the disk and a linear 
longitudinal potential distribution on the convex 
surface, the equations and tabulated coefficients 
given in this paper permit rapid hand computation 
of thermal conductivity values from experimental 
data. If boundary conditions other than those of 
eq (9) are required or indicated, thermal conductivity 
values should be computed from (12), rather than 
(13), with the appropriate coefficients being given 



5 Cf. equation (12), page 514, reference [3]; in the present paper this equation 
has been reduced, by means of recurrence relations, to a simpler form than that 
given originally by lloch et al. 
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by (6), (7), and (8). Equation (12) is completely 
rigorous for the axially-symmetric boundary con- 
ditions (4a), (4b), and (4c). Proof of adequate 
compliance with these boundary conditions would 
appear to be indicated for any given test method. 
Hoeh et al., [3, 9] have discussed the influence of one 
particular type of boundary condition, other than 
linear, on the convex surface of the disk. Other 
boundary conditions can be handled by evaluation 
of the appropriate integrals given in this paper. 



The author thanks B. A. Peavy for checking the 
mathematical analysis given in this paper. 
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